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Abstract 

A theoretical model is proposed for the (0,0,1) superlattice manganite system 
(LaMn03)n(SrMn03)m- The model includes the electron-electron, electron-phonon, and co- 
operative Jahn- Teller interactions. It is solved using a version of single-site the dynamical mean 
field approximation generalized to incorporate the cooperative Jahn- Teller effect. The phase 
diagram and conductivities are calculated. The behavior of the superlattice is found to a good 
approximation to be an average over the density-dependent properties of individual layers, with 
the density of each layer fixed by electrostatics. 

PACS numbers: 68.65.Cd, 71.10-w,73.20-r,73.40.-c 
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I. INTRODUCTION 



An exciting recent development in materials science is the ability to fabricate atomically 
precise multilayer structures involving transition metal oxides [l, 2], s], ^,0,0]. Multilayer 
structure based on simple semiconductors such as Si and GaAs/AlGaAs give rise to wide 
range of striking physical phenomena such as the integer and fractional quantum Hall effects 
0, Q] and are important for many classes of devices. It is therefore natural to expect that 
transition metal oxides, which display a far richer variety of phenomena than do semicon- 
ductors jsl, will in heterostructure form yield an even more diverse set of new phenomena. 



The colossal magnetoresistance 



studied intensively 
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CMR) rare earth manganites Lai-^Sr^-MnOs have been 



14, 



15l | and provide an important model system. These 



compounds display a wide range of phases with characteristic orbital, magnetic and transport 
signatures 10|, lUl, UM, llTI], providing opportunities both for creating many effects in the 



superlattice context and for detecting phases that may be created. Mariganite superlattices 



have now been fabricated and are being studied experimentally 0, 0, 16| but have not yet 
received much theoretical attention. In this paper we present a theoretical analysis of these 
interesting systems. Our findings may be summarized by the following 4 rules: (1) the 
layer charge distribution mainly depends on the electrostatic interaction and has very weak 
dependence on the temperatures and on whether or not the system has orbital or magnetic 
order; (2) once the layer charge distribution is given, the orbital order at each layer is 
essentially determined by the bulk behavior at the same density: the propagation length 
of the orbital order along the superlattice is less than a lattice constant; (3) the inter- 
layer magnetic coupling is ferromagnetic (FM) except between layers with densities close 
to 1, where it is anti-ferromagnetic (AFM); (4) the in-plane conductivity is essentially the 
average of bulk conductivities for the given layer charge distribution. These results indicate 
that the main effect of superlattice structure is simply to produce a layer-dependent charge 
distribution. Once the layer densities are known, each layer behaves basically according to 
the bulk phase diagram for that density with only weak proximity effects. 

The rest of this paper is organized as follows. First we define the superlattice and discuss 
the interactions considered in the bulk model, then we present and discuss results for a 
particular superlattice consisting of four layers of LaMnOa and one of SrMnOs in detail. 
From these and related calculations we infer the rules outlined above. Finally we use the 
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rules to discuss recent experiments j^, l5| arguing that the data still may not be displaying 
intrinsic behavior. We propose experiments on a variant of the superlattice which may be 
more definitive. 

II. DEFINITIONS, HAMILTONIAN, METHODS, AND PHASES 
A. Definition of Superlattice 
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FIG. 1: Upper panel: Schematic representation of a (LaMn03)4 (SrMnOa)! superlattice with Mn 
manganese and O oxygen ions denoted by shaded and open circles respectively. Lanthanum and 
Strontium ions labeled by La/Sr reside in the center of cubes defined by Mn. Lower panel: Upper 
part - Projection of (4,1) superlattice onto x — z plane, with the three symmetry inequivalent Mn 
sites labeled A,B,C. Middle part - a numerical representation of charge density and amplitude 
of orbital order computed for screening parameter a = 0.3 at temperature indicated. Lower 
part - pictorial representation of magnetic order (order parameter direction indicated by arrows) 
calculated using screening parameter a = 0.3 for two temperatures as indicated. 

The high temperature structure of LaMnOs is a slightly distorted version of the ABO3 
perovskite. It may be thought of as a simple cubic lattice with Mn on the cube vertices and 
La in the body centers. The superlattices of interest here are composed of the same cubic 
lattice of Mn, with a superstructure defined by a periodic replacement of La by Sr. To specify 
the superlattice we must specify the direction along which the La/Sr sites alternate and 
the way the La and Sr ions are arranged. We will consider a (0,0,1) (LaMn03)m(SrMn03),i 
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superlattice, abbreviated as LmSn or {m,n), defined by m La planes perpendicular to the 
(0,0,1) cube axis, followed by n Sr planes. The entire structure has a periodicity (m + n) in 
the z (coordinate) direction while maintaining the original lattice translational symmetry in 
X and y. The upper panel of Fig([T]) shows a L4S1, (4,1) superlattice. We shall refer to each 
X — y plane as a layer. 



B. Hamiltonian and Parameters 



The degrees of freedom needed to describe the physics of manganites are [Uj the two 
Cg-like orbitals which make up Mn conduction band, an electrically inert classical core spin 
representing electrons occupying orbitals at Mn sites, and the three even-parity MnOe 
distortion modes Qo, Qx, Qz- The Hamiltonian for the manganite superlattice is the Hamil- 
tonian for bulk manganites, supplemented by Coulombic terms representing the potential 
arising from the pattern of the La and Sr ions, thus H = Hsuik + Hcoui- The bulk Hamil- 
tonian is HbuIU = Hband + HeE + Huund + Hel-lat + Hiattice where 

Hband = ^ ^k,ab,a^\a,a%b,a (l) 
k,ab,a 

with e^^^^ = —t{eQe + SzTz + SxTx)ab where e is the unit matrix, f Pauh matrices and 
Eq = cos kx + cos ky + cos k^, Sz = cos k^ — |(cos kx + cos ky), and Sx = ^(cos kx — cosL, . 



a, b label orbitals, i,j sites, and a spins. From our previous band structure calculation 18| . 
the hopping parameter t is 0.65eV which defines the energy unit. 

Hee = {U - J)nin2 + ^ '^iA'^hl + J{c\,^c\,f'i,l'^'iA + ^■^■) - 2J5i ■ ^2 (2) 



j=l,2 



with 5i(2) = '?a/3c|{2),a'^l(2: 

priate for the manganite 



19|. 



The values U ~ 2.5eV and J ~ 0.5eV are found to be appro- 



HHund — —Jh'^^ Si ■ c\^^aal3Ci,^ (3) 



with Jh ~ 1.5eV and 15*1 = 1. An additional antiferromagnetic core-spin core-spin inter- 
action exists in the material. This interaction is crucial for x > 0.5 but affects only minor 
details of the x < 0.5 regime of interest here, i.e. it changes the phase boundaries but does 
not lose any exhibited phases. Including this interaction involves too great a computational 
expense, so we do not include it here. 



For the electron-lattice interactions [12] we only include the Jahn- Teller coupling 



(4) 



i,a,b 



with Qi 
coupling 



Qi,X} Qi 



Qi{sm29i,cos29i) {0 < 9 < vr). We neglect the breathing mode 



12 1 in the current calculation. 



For the lattice Hamiltonian we include an elastic energy term from adjacent Mn-0 and 
Mn-Mn force constants 



lattice 



15l . |20| and the cubic term noted by Kanamori [21 1 
E ^bQ^,aQj,b - A J2{3Ql - QlQ.,.) 



2K" 

ij,ab *J 



(5) 



with i,j labeling site while a,b distortion modes. 

ectrons are concerned, La and Sr ions act as +1 and neutral point charges 



As far as e 



respectively 
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23|. In the superlattice, the distribution of those cations is patterned, and 



the Hamiltonian from electrostatics is 

Hcoul = 



1 e^riirij 1 
^ + 



2e m 



2e |i?f« - 



e\fi - Rf^l 



(6) 



with rii = 'Yl,ua^ta(T^i,a,cT the occupation number at Mn site i. ri and -Rf" are positions of 
Mn and La in ith unit cell, and e is the dielectric constant of the material. To describe the 
magnitude of this interaction, we define the dimensionless screening parameter a = e'^/aet 
which controls the charge density distribution. The value of a has not been determined, but 
its order of magnitude may be estimated from the hopping parameter t ~ 0.65eV, lattice 



constant a = 4A and typical value of dielectric constant e ~ 10 24 1 to be a = 0.2. 



C. Method and Results from Bulk Model 



We now briefly discuss the method of solving the rnodel. To solve the model we use 
;he single-site dynamical mean field theory (DMFT) ^] with semiclassical impurity solver 
26 1 . The intersite lattice coupling (cooperative Jahn- Teller effect 0, Q) is taken into 



account within the single-site DMFT by the mean field approximation, i.e. the local impurity 
problem contains an effective interaction from distortions at neighboring sites 15|. For 
the term introduced by the superlattice Hcoui, we adopt the Hartree approximation whose 



overall effect is to generate a self-consistent potential at each Mn site 
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23| . The optical 



conductivities are computed in the same manner as described in Ref[lJ], and the DC values 



are obtained by taking the 
in Ref|l5|. 



hmit. Sepcifics of the DMFT computation are presented 



Calculated Phase Diagram 
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FIG. 2: The calculated bulk phase diagrams. To make comparison between our calculations and 
experiments, the calculation done at T = O.lt correspond to roughly 250K in the experiments [l^. 



The first step to study the superlattice is to understand the parental bulk system. Fig([2]) 
shows the theoretically calculated bulk phase diagram. Fig(l2]) is diff'erent from what we 
presented in Refjl^ in phase boundaries because here we neglect the AF magnetic coupling 
between Mn core spins. This simplification does not lose any phases for the doping range 
X < 0.5 but does lose phases for x > 0.5, therefore we only present bulk results for x < 0.5 
and only consider {n{> 2),1) superlattices where all layers have Cg electron occupations 
larger than 0.5. Our bulk calculation overestimates the transition temperatures but gives a 
reasonable description of the ordering of phases. For this reason we focus on the exhibited 
phases in superlattice but not the transition temperatures on the numerical values. 



D. Planar Long-ranged Orders 



Similar to earlier findings [22|, |23| , the superlattice results are most naturally presented 
by first describing the charge density and exhibited order(s) at each x — y plane, and the 
superlattice phase are therefore the planer phases plus the relations between them. To 
facilitate later discussion, we summarize here the three planar long-ranged orders found in 
our calculation: a (vr, vr) staggered Qx, a uniform Qz, and a ferromagnetic (FM) order. The 
first two are orbital orders; the ordering pattern are shown in FigQ. For the staggered 
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FIG. 3: The orbital orders, (a) The staggered Qx order and (b) the uniform Qz order. For the 
staggered Qx order, all octahedra in the plane undergo Qx distortions with the long-axis (double 
arrows) of each distorted octahedron differs by 90 degree from its adjacent sites. For the +{-)Qz 
order, all oxygens are in the middle of neighboring Mn, but the Mn-0 distance in z direction is 
longer (shorter) than that in x and y. 

order, all octahedra in the plane undergo Qx distortions with the long-axis of each distorted 
octahedron differs by 90 degree from its adjacent ones, as shown in Fig([3])(a). For the +(- 
)Qz order, all oxygens are in the middle of neighboring Mn, but the Mn-0 distance in z 
direction is longer(shorter) than that in x and ?/, as Fig([3])(b). 

An equivalent way to describe the local orbital order is to use the ground state \6) = 



cos^|32;^ — r^) + sin^|x^ — y"^) 151] of the Jahn- Teller coupling (Eq(j4])) with Qi given by the 
distortions. In this language the (vr, tt) Qx order corresponds to an alternating \6 = vr/4), 
\6 = — 7r/4) orbital configuration and the uniform +{—)Qz a uniform \6 = 0) {\6 = 7r/2)). 



III. RESULTS FOR (4, 1) SUPERLATTICE 

In this section we present our calculated results for the (4, 1) superlattice, a representative 
case that captures most of the relevant phenomena. We consider two different values of a; 
a = 0.09 corresponding to relatively delocalized electrons and a = 0.3 corresponding to 
tightly confined electrons. 

For L4S1 each supercell has 3 symmetry-inequivalent layers, labeled as A B C where A is 
sandwiched by La, B by La and Sr, and C by Sr, as illustrated in Fig([T])(b). Fig(j4]) presents 
the L4S1 layer charge distribution for a = 0.09 (corresponding to a relatively delocalized 
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FIG. 4: L4S1 charge distribution for a = 0.09 at 0.16t (above any ordering temperatures), 0.12t 
(orbital ordering but not magnetic ordering in super lattices), and 0.04i (both orbital and magnetic 
order presented). 



electrons) at temperatures 0.16t (above any ordered temperatures), 0.12t (orbital ordering 
but not magnetic ordering presented), and 0.04t (both orbital and magnetic order). The 
variation in charge distribution for different temperatures is smaller than 5%. 

The staggered order. Table I lists the layer charge density, the computed SL order, 
and the bulk order at the given layer charge density of the (4,1) SL for a = 0.3 and 0.09 
at T = 0.12t. We see that the SL calculation is very close to the bulk (< 1%) which means 
that the Qx order is mainly determined by its layer charge density and the proximity effect 
is very weak. 





a=0.3 


a=0.09 


site label 


A 


B 


C 


A 


B 


C 


Charge Density 


0.99 


0.99 


0.52 


0.92 


0.90 


0.64 


ira|SL[t] 


1.17 


1.17 





1.01 


1.0 





|(Q^)|Bulk [t] 


1.18 


1.18 





1.02 


1.0 





|(Q^)| SL [t] 


-0.28 





-0.14 


-0.18 


-0.03 


-0.08 


Bulk [t] 


-0.32 


-0.32 





-0.22 


-0.2 






Table I: The amplitudes of the staggered Qx i\{Ql)\), uniform (KQ^)!) 
orders for the L4S1 superlattice and for the corresponding bulk values 
at same densities for T = 0.12t 



The very weak orbital order proximity effect for the staggered Qx order has a straightforward 
physical origin. Rotation symmetry about the bond in the z direction means that neither 
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nearest neighbor (NN) hopping nor nearest neighbor atom-atom forces can distinguish +Qx 
order from —Qx order. One consequence is that within this approximation (tt, tt, tt) 
and (vr, vr, 0) Qx order are degenerate implying that the in-plane (vr, vr) Qx order does not 
have any dispersion in z direction and thus no inter-layer coupling (no proximity effect) for 
this order. Second neighbor hopping [15] or certain classes of shear elastic forces [28.] do 
distinguish the two but the differences are small (of the order of meV). We will give a more 
detailed discussion in the discussion section. 

The uniform Q^ order. The 3rd and 4th rows in Table I list the computed SL Q^ order of 
the (4,1) SL for a = 0.3 and 0.09 at T = 0.12t and the corresponding bulk results at the same 
density. These results can be understood as follows. There are two apparent sources affecting 
the uniform Q^ order - the anharmonic term in Hiatuce and the SL effect. The net SL effect is 
to produce an electrostatic like force, i.e. oxygens are pushed toward to Mn layer with less Cg 
density. One bears in mind that the origin of this force is because the octahedral distortion 
changes the Mn-0 hybridization thus affects the population distribution of eg orbitals, and 
has little to do with the actual static Coulomb energy [l^. For the parameters considered 
above. Layer A is sandwiched by Layer B having similar charge densities (1 in this case), 
therefore the SL effect is weak for A and the main contribution is from the anharmonic 
effect. Indeed the Qz order of Layer A is -0.28t, very close to the bulk result -0.32/!:. For 
layer B, besides the bulk effect which causes a -Qz order, the charge inhomogeneity results 
in a +Qz order since oxygens are pushed away from layer B and the in this case these two 
effects nearly cancel each other. For layer C, only the SL effect contributes and the result 
is a 0.14t ~Qz order. We see that at interfaces the bulk and SL effects are of comparable 
strengths therefore the uniform Qz order has significant proximity effect. 

The magnetic order. Now we discuss the magnetic order for a = 0.3 and 0.09. For 
a = 0.3, only layer C has FM order at T = O.lt. When lowering the temperature to 0.04t, 
all layers are 2D FM. The interlayer coupling between A and B (with layer densities roughly 
one) is AF, and all others are FM as shown in Fig([T])(b). For a = 0.09, one obtains a 
smoother charge distribution and at low temperature all inter-layer magnetic couplings are 
FM. Those results suggest that the interlayer magnetic coupling is generally ferromagnetic 
except that between layers with densities close to one. The FM coupling between layers 
can be understood from double exchange mechanism where FM arrangement maximizes the 
kinetic energy and is energy favored when the conduction band is not full. The AF coupling 
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between = 1 layers has the same origin as the bulk ISj where the system gains gap energy 
by arranging core spins antiferromagnetically. One interesting consequence of these results 
is that if the charge density is sufficiently sharp enough to introduce AF inter-layer coupling 
{a = 0.3 for example), the superstructures with odd number of La have an odd number of 
AF bonds so the magnetic lattice has twice the period of the structured one. For example 
for the L3S1 superlattice, we expect the period including the magnetic order is 8 and adding 
up magnetic orders from all layer leads to zero total magnetization. 



IV. SPECTRAL FUNCTIONS AND CONDUCTIVITIES FOR DIFFERENT SU- 
PERLATTICES 

In this section we shall first establish the rule implied from our superlattice calculation 
for the in-plane conductivity and also discuss the out-of-plane one. Based on the established 
rule we discuss the temperature and superlattice dependence on the DC measurements. 



A. The rule for the conductivity 

We found that once the layer charge densities rii are specified, the excitation spectrum 
at each layer behaves very close to the bulk at the same density. We use the simplest 
superlattice L2S1 to demonstrate this point. For (2,1) SL there are only two symmetry- 
inequivalent sites, one between La labeled as A while the other between La and Sr B, as 
shown in Fig(l5]) (a). For a = 0.3 and T = 0.12t we found ~ 1 and ub ~ 0.5. Fig(l5]) 
(a) /(b) shows both the calculated SL spectral functions for layer B(A) and the bulk for 
A^ = 0.5(A^ = 1) at the same temperature. 

Since the in-plane conductivity af^ is dominated by in-plane hopping processes (moving 
electrons from one site to the other within the same layer) and the spectral function for 
layer with density rii is very close to the bulk at the same density, crf^ should be very close 
to the bulk-averaged value o"^^''^ which is define as 

N 

^'T^Y.^lTin^)/N (7) 

i=l 

To see how accurate this rule is, we list in Table III the calculated SL DC conductivity 
crf^ and the bulk averaged a^^''^ according to the obtained charge distribution at T = 0.1 2t 
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L2S1: 




• — « Site A, majority oibiliil 
Q — o Site A, minority orbital 

N=l, Bulk, majority orbital 

N=l, Bulk, minority orbital 



Frequency (t) 




(a) (b) 

FIG. 5: Spectral functions for the (2,1) superlattice compared to bulk spectra at density equal to 
layer density for a = 0.3 and at T = 0.12t. (a) The superlattice spectral functions for layer B 
and bulk at = 0.5. (b) SL spectral functions for layer A and bulk at = 1. To facilitate the 
comparison, the bulk results are shifted OAt in frequency. The peak-peak distances are essentially 
the same for bulk and SL calculations. 

where the unit for axx is e^/{ha) ~ 6 x 10^(i7cm)~^ with a the lattice constant (~ 4A). 



a 


0.3 


0.09 


0.3 


0.09 


0.3 


0.09 


-^bulk 


0.083 


0.115 


0.064 


0.093 


0.58 


0.082 


^SL 

XX 


0.077 


0.105 


0.067 


0.093 


0.6 


0.081 


SL 


L2S1 


L2S1 


L3S1 


L3S1 


L4S1 


L4S1 



Table II afj" and o"^"''^ for different superlattices 

As a concrete example we compute bulk-averaged o"^^'^ for L4S1, a = 0.09, T = 0.12t in 
details. From FigQ we see that there are two layers with density 0.9 (layer B), two layers 
0.64 (layer C), and one layer 0.91 (layer A) for those parameters. Now we consult the bulk 
results for DC conductivity at T = 0.12t and find that a^J'J'' at =0.9, 0.64, 0.91 are 
0.055, 0.12, 0.05 respectively. The bulk averaged is thus (2a^^''=(n = 0.9) + 2af,^,'^(n = 



0.64) + a 



bulk I 



n 



0.91))/5 roughly 0.081. As listed in Table II, the difference between cr^"''' 



and a^^ is roughly 10% for L2S1, and is smaller (~ 2%) for L3S1 and L4S1. FigQ compares 
the in-plane optical conductivities with the bulk-averaged for the L4S1 at a = 0.3, T = 0.12. 
We see they match quite well (< 10%) (both peak positions and peak amplitudes) for all 
frequencies. 
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FIG. 6: Calculated SL in-plane optical conductivity and the corresponding bulk-averaged one for 
the L4S1 at a = 0.3, T = 0.12t. As a reference we also plot the bulk axx{^) foi" N = 1 and = 0.5 
(curves without ligands) at the same temperature which are used for computing the bulk-averaged 
conductivities. 

Unlike the in-plane case, the out-of-plane conductivity a^^ involves inter- layer hopping 
processes (removing electrons from layer i and adding them to layer i ± z). Since different 
layers experience different static Coulomb potential, there is no simple relation between af/' 
and {rii}. Our previous study 22| of a simpler model system, namely the one-orbital double- 
exchange superlattice, shows that there are peaks in (J^^{uj) directly corresponding to the 
potential difference between layers from which the a is determined straightforwardly. How- 
ever for the current problem the potential difference is mainly to produce layer- dependent 
local spectra and are not directly related to peaks in erf/'. Nonetheless af/' is more sensitive 
to the screening parameter than erf/ and can be used to constrain the values of a. We will 
return to this point in the discussion section. 



B. Temperature dependence of DC resistivity 

Now we discuss the temperature dependence on the DC resistivity p. Fig([7j) shows p(T) 
for a = 0.3. We observe that around T = O.llt dp/dT changes sign implying the superlat- 
tice goes from high T insulating phase to low T metallic. This downturn in p(T) coincides 
with the interface Curie temperatures which are also marked Fig(I7j). Across this temper- 
ature the interfaces go from PM/Bad Metal to FM/Meta^ 
metallic behavior. From double-exchange mechanism jl4 
always accompanied with a insulator/metal (bad metal/metal at least) transition, the sign 
change in dp/dT as a function of temperature is quite general for all values of a, but is more 



accounting for the superlattice 
Q, the PM/FM transition is 
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0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 
Temperature (t) 



FIG. 7: The DC resistivity as a function of temperature for a = 0.3. X indicates the calculated 
interface Curie temperature for each superlattice. The vertical dashed line marks a rough tem- 
perature where dp/dT changes sign. As a reference we also present the bulk result for N = 2/3 
{x = 1/3). 

pronounced for large a. 

C. Superlattice dependence of DC conductivity 

Since for small a, the superlattice behaves just like the bulk material, we only focus on 
the a which leads to the sharp charge distribution, i.e. a = 0.3. Our result implies that the 
SL conductivity is proportional to the interface density at low temperature (the interface 
density is defined as the ratio between the number of Mn sandwiched by La and Sr and that 
of total Mn layers, for L2S1 the interface density is 2/3, for L4S1 5/2). This statement is 
a direct consequence of the rule crfj" ~ a^^''^: for sharp charge distribution, ~ 1 layers 
are insulating at low T and only interface layers (A^ ~ 0.5) are conducting. To demonstrate 
this statement from the SL calculation, we define r(T) = cr^^'^^ (T) /cr^J'^^ (T) and compare 
it with the ratio of interface densities between L4S1 and L2S1 which is 0.6. The following 
table presents r{T) for several temperatures: 



T/t 


0.15 


0.1 


0.07 


0.04 


r(T)"^'^''^ 


0.96 


0.65 


0.55 


0.57 



We see that r{T) is indeed very close to 0.6 below the interface Curie temperature (~ O.llt). 
Finally we emphasize that from our calculation, the statement "SL conductivities are only 
from interfaces" is true only below temperature and with sharp charge distribution. 
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V. DISCUSSION 



In this section we first summarize our results by providing rules deduced from our calcu- 
lations, then discuss how to determine the layer charge distribution and thus the screening 
parameter experimentally. Finally we give a more detailed analysis on the proximity effect 
of the orbital orders. 



A. Rules implied by the calculation 

The exhibited phases at each layers are the combined effects of bulk Hamiltonian and the 
charge inhomogeneity induced by the SL. Here we summarize our calculations by stating 
the following rules which governs the displayed phases at each layer. 

l.The charge distribution is mainly determined by the screening parameter a, and is not 



sensitive to the temperatures and orders 



22 



23|. 



2. The staggered Qx order at each layer follows the bulk behavior at the corresponding 
layer charge density density. 

3. The uniform order is caused by two sources. First (bulk effect), a uniform —Qz 
order is induced with the presence of the staggered Qx order. Second (SL effect), the charge 



inhomogeneity induces an electrostatic like force pushing oxygens closer 



;o the layer with 



lower density. Note that the origin of this force is the Mn-0 antibonding 

4. The magnetic order. The inter-layer coupling is FM except that between layers with 
densities close to 1. 

5. The in-plane conductivity of the SL is essentially the average of the bulk conductivities, 
i.e. for a given charge distribution {ni,n2, ..n]y}, a^^ = a^^^^ = J^iLi '^l^'''' {''^i) / ^ while we 
do not have simple rule for out-of-plane conductivities. 



B. Determining the charge distribution and screening parameter 

Our results indicate that the key quantity for the SL is the charge distribution and 
therefore determining the screening parameter a is very important. In principle the charge 
distribution can be determined by measuring the change in manganese valence in the super- 
lattice but such measurements have not to our knowledge been performed on manganite 
superlattices. Here we propose two measurements to constrain the values of a. The first 
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FIG. 8: In-plane (a) and out-of-plane (b) optical conductivities for L2S1 for a = 0.09 (with ligands) 
and 0.3 (without ligands) at T = 0.12t. 

proposal is to grow superlattices with odd number of LaMnOa layers (L3S1 for example) and 
to measure the total magnetization at low temperature. If the charge distribution is sharp, 
then the ferromagnetic interface layers are antiferromagnetically coupled and the layer has 
zero net magnetization, whereas if the charge distribution is broad, the system will simply 
be ferromagnetic. This method suffers from a potential experimental disadvantage: lack of 
perfect oxygen stoichiometry could produce ferromagnetism in the La-rich regions even in 



the presence of a sharp charge distribution 4. |29|. 



The second one is to measure the optical conductivity. The basic idea is that one should 
be able to decompose the in-plane optical conductivities into sum of bulk results. This 
measurement is free from the interferences of the impurity scattering or extra oxygens in the 
superlattice. Fig([8])(a) shows the in-plane L2S1 optical conductivities cr^x^^i^) = 0-3 

and a = 0.09 at T = 0.12t (~300K and with interfaces being slightly ferromagnetic). We 
see that for a = 0.3, the cr^x^^iuj) is not peaked at zero frequency because 1/3 of the 
contribution is from Mn layers with density close to one. On the other hand for a = 0.09, 
'^xx^^i^) behaves more like a bulk material. 

The out-of-plane conductivities contain information of the potential differences at dif- 
ferent layers. As discussed in section IV. A, as for the azz is concerned, the main effect of 
the layer-dependent potentials is to produce layer-dependent local spectral functions which 
makes different from a^x at high T non-ordered phase. Therefore the anisotropy in a 
above ordering temperatures is a very straightforward estimator of how confined/delocalized 
the electrons are. Fig(IH])(b) shows the out-of-plane conductivity crff'^^(u;) for the (2,1) su- 
perlattice. As expected larger a results in larger anisotropy in cr. For the numerics, because 
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the potential difference is not directly related to the peak in the out-of-plane optical conduc- 
tivity as in the model studied in Ref 22], one has to compare the experimental data with 
theoretical results. 



C. Orbital order proximity effects 

Our calculation indicates that orbital order associated with the Qz mode (octahedral 
distortion with long bond pointing along z direction) is reasonably efficiently transmitted 
from one layer to the next, whereas orbital order associated with (octahedral distortion 
with long axis in x-y plane) is not transmitted at all. The approximations made in this 
paper amount to retaining only nearest neighbor interactions in both the electronic and 
lattice sectors; in this approximation rotational invariance about the bond in the z direction 
ensures that Qx order is not transmitted. In this section we estimate the extent to which 
terms not included in our model may change this result. 

The transmission of Qx orbital order implies an energy difference between states in which 
Mn ions in adjacent planes have the same or opposite amplitudes for Qx order. To estimate 
this we consider the energy difference between (tt, tt, 0) and (tt, tt, tt) orbital order. There 
are two contributions to this energy: electronic and elastic. 

We begin with the electronic term. A given lattice distortion selects one locally favored 
state of each site i, which we refer to as 6'f >= cos6'j|3z^ — > +sin6i\x^ — > and 
also the orthogonal excited state Q\ >. For the +Qx distortion the locally favored state is 
\6q^ >= ^|3z^— > +--i=|x^— > and the locally disfavored state is 6^^^ >= -^\3z'^—r^ > 
— — >; for the —Qx distortion the roles of 6g and 6'^ are reversed. 

The (tt, tt, 0) Qx and (tt, tt, tt) Qx orders produce exactly the same local distortions; an 
electronic contribution to the energy difference between these two orders can therefore only 
arise from a difference in the electronic hybridization energies for the two states. We now 
argue that if only nearest-neighbor hopping is present, there can be no hybridization energy 
difference between the two states. The in-plane hopping is trivially the same for the two 
states; a difference can therefore only arise from the z direction hopping. The hybridization 
is given by a 2 x 2 matrix; in the nearest neighbor hopping approximation the matrix is 
degenerate in all three directions with only one non-zero eigenvalue, t. For hopping in the 
z direction the only non- vanishing matrix element is the one connecting \3z^ — > orbitals 
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on the two sites implying that the matrix elements between any combination of \0q^ > 
and|6'Q^ > is t/2, so that the nearest neighbor hopping does not distinguish between the 
two states. Including the 2nd NN couplings can lift this degeneracy. The energy difference 
can be estimated by comparing the energy gain caused by the virtual processes involving 2nd 
NN hoppings (super-exchange) and is done in Ref The LDA calculation [isl suggests 
the 2nd NN hopping is roughly 0.035eV, leading to an energy difference of the order of meV 
for these two orders. 

We next consider contributions arising from elastic forces. In a ball and spring model, 
we have found that oxygen-oxygen or three-body forces are required to propagate Qx order. 
For these situations, spring constant models are not reliable and band theory calculations 



of phonon stiffness are required. Available information on these forces is presented in |28j : 
the result is that the energy difference is of order meV, roughly the same magnitudes as 
the electronic contribution. With these estimates we conclude that our model captures the 
main physics of manganites and the proximity effect for the planar staggered Qx order is 
indeed very weak in reality. 

VI. CONNECTIONS TO EXPERIMENTS 

The temperature dependence of the DC resistivity (p) and magnetization (M) of 



(LaMn03)2n(SrMn03)n (the (2n,n) superlattice in our notation) have been measured 
In this section we discuss the connection between our calculation and these measurements. 



A. DC Resistivity 

For the resistivity measurement, as the temperature is decreased from 400K, the resistivity 
p is found to increase until a superlattice-dependent temperature is reached. This 
"downturn" in p(T) is observed for all {2n, n) superlattices. The temperature at which the 
downturn occurs (T^) is n dependent (larger n lower Td) and roughly 200-350i^ (Fig(2) 
(a)-(c) in |4|], Fig(2) in [5|). Above T^, the magnitude of p has a very strong dependence on 
n. When the temperature is further reduced below lOOi^', the resistivity increases again for 
n > 3 superlattice while keeps on decreasing for n < 3. 

Calculation finds two typical behaviors corresponding to weakly and strongly bound 
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n ( Superlattice (2n,n) ) 



FIG. 9: The in-plane DC resistivities from experiment at T = 150K [5] (solid line) and from Rule 
5 at T = O.OSt with sharp charge distribution. 

charge distributions For both cases we expect a downturn in resistivity which is associ- 
ated with the interface FM/PM transition with the downturn temperature essentially n 
independent. At very low temperature our theory predict metallic {dp/dT > 0) behavior. 
The n dependence on the magnitude of p are different for these two limits: for the weakly 
bound charge distribution p has very weak n dependence while for the strongly bound p is 
proportional to the interface density and thus oc n. 

The results according to our rules are in reasonable agreement with the n < 3 superlattices 
because they remain metallic in the low temperature. For long period superlattice except 
the existence of the changing sign in dp/dT, our theory is not compatible with experiments 
at quantitative level because of the exhibited n-dependence on p and T^, and because of the 
low T insulating phase. Essentially the experimental measured p depends too strongly on n 
compared to our theoretical prediction as shown in FigQ. We think these inconsistencies 
are caused by the ra-dependent interface quality. 



B. Magnetization 

For the magnetization measurement, Ref ^ (Fig(2) (d)-(f )) shows that all superlattices 
exhibit a net magnetization below some onset temperature T/ (the superscript / denotes 
interface). T/ depends on n (decrease in T/ from 300K to 200 when increasing n from 2 to 16 
) and is very close to the downturn temperature in resistivity T^. For the n = 2 superlattice, 
the magnetization increases and saturates when lowering the temperature. For n = 4 n = 16 
superlattices, the magnetization increases slowly right below T/, then more rapidly around 
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150K. When the system is further cooled without applying a magnetic field (zero field 
cooling, ZFC), the magnetization reaches a maximum (roughly the same value ~ 0.4/iB/Mn 
for both n = 4 and n = 16) around llOK, before eventually decays to zero. If the system 
is cooled in an applied magnetic field (field cooling, FC), the n = 2 superlattice behavior 
is not changed but for both n = 4 and n = 16 the magnetization saturates at 0.6/iB/Mn 
approaching zero temperature. We believe that the FC measurements reveal more clearly 
the intrinsic behavior and will therefore focus on the FC data but will comment on the ZFC 
data in the end. 

Our calculation predicts that ferromagnetism is associated with the interfaces. There 
is a critical thickness which depends on the width of the charge distribution (value of a). 
For systems with n greater than this critical value (which may be for very sharp charge 
distributions), the magnetization resides only on the interfaces, so one expects an onset at 
the interface (comparable to the bulk value) and a saturation value scaling with the 
interface density (with coefficient depending on a). 

Now we discuss the magnetization (the FC data) based on our calculation. First our 
theory agrees reasonably well with the n = 2 (L4S2) superlattice assuming the critical 
thickness is > 3 so the LaMnOs layers are ferromagnetically ordered (the n = 2 superlattice 
saturates at 2.8/iB/Mn at low temperature indicating LaMnOs layers are ferromagnetic). 
However, the behavior of the n = 4 and = 16 superlattices is not consistent with our 
theory, except in the qualitative sense that that resistivity downturn and magnetic onset 
temperatures coincide. The strong n-dependence of the magnetism-onset temperature T/ is 
not found in the calculation, where ferromagnetism is an interface phenomenon and therefore 
not strongly n-dependent. Further, the observed rapid growth in M below 140i^ and the 
saturation at about 1 /i^ displayed in the long period {n > 6) superlattices is not found 
in our theory, which instead implies antiferromagnetism below 140i^ and a saturation 
magnetization proportional to the interface density. 

We believe that the inconsistency is a consequence of oxygen stoichiometry that extra 
ox yge ns reside in La-riched regions so that LaMnOs is effectively a slightly doped manganites 
ji, 29|. Based on this assumption and our rules, we now explain the FC data. Our rule states 



that each layer behaves as a bulk manganite at the given charge density. 



close to but not equal to one, the system displays canted A-AF order 



or the density 



nj around 140K 



which has a small but non-zero net magnetic moment. Within this picture, the constant 
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but small saturation magnetization for n > 6 superlattices is the net magnetization from 
the canted A-AF order, and the rapid growth around 140K which is caused by the forming 
of the canted A-AF order in the La-riched regions. The remaining issue is the n dependence 
on the interface Curie temperature. We believe this is caused by the n-dependent interface 
quality. Since the main energy gain for PM/FM transition is the kinetic energy (double- 
exchange), we expect the bad interface quality tending to localize electrons will reduce the 
Curie temperature. This is qualitatively consistent with the observation. 

Finally we comment on the zero field cooled data. The most striking feature for ZFC data 
is that for the ri > 4 superlattice the magnetization starts to decrease around llOK upon 
cooling. We associate this reduction in M with the G-AF order developed in the SrMnOa 
layers. The forming of G-AF order in SrMnOa weakens or temporarily kills the ferromagnetic 
coupling between ferromagnetic interfaces therefore reduces the total magnetization. Note 
that the Neel temperature for CaMnOa is llOK 16|. 



C. Interface Qualities and Proposed Experiments 

Although the qualitative interpretations can be made, our rules are not compatible with 
experimental data at quantitative level, in particular the resistivity measurement. We be- 
lieve this discrepancy is mainly caused by the n-dependent interface quality which has two 
possible causes. The first one is the mismatch in lattice constant. Since the LaMnOs and 
SrMnOs have slightly different lattice constants, the interface either tolerates very large 
strain forces or more probably creates some defects to compensate the mismatch. The 
presence of defects increases the scattering sources and thus the resistivity which is not con- 
sidered in our model. The second one is the magnetic frustration. Since LaMnOa, SrMnOs 
at low temperature display A-type and G-type AF orders respectively, interfaces undergo 
the magnetic frustration which suppresses the tendency of FM order and again increases the 
resistivity. 

To reduce these extra variables in experiments, we propose one should conduct experi- 
ments on (n, 1) SL. For this (n, 1) series SL, since no Mn are sandwiched by Sr layers we 
expect the effects due to the lattice constant mismatch should be relatively small. Perhaps 
more importantly complications from the magnetic frustration caused by SrMnOs AFM 
regions will not enter the problem. Based on these arguments the n dependent interface 
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quality is expected to be strongly reduced and the direct comparison between experiments 
and theory becomes easier. 

VII. CONCLUSION 

We have studied the (LaMn03)„(SrMn03)i (n, 1) superlattice. We found that the charge 
distribution, determined by the screening parameter a, is the key quantity and propose that 
the optical conductivity measurement of L2S1 can fix it. Once the charge density at each 
layer is known, the phases essentially follows the bulk at the corresponding density. General 
rules for phases are given. We propose the measurements on low temperature magnetization 
of the L3S1 SL and on the optical conductivity help constrain the range of a. For a — 
0.3 which results in sharp charge distribution, we predict at low temperature the (n, 1) 
superlattice has no net magnetization with period 2 x (n+ 1) for odd n. Finally we comment 
on recent experiments based on our calculations and infer that the main discrepancy is caused 
by the sample-dependent interface quality. We suggest the (n, 1) superlattice is actually a 
physically simpler and cleaner system to study. 
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